



use "stata data/regsave/A1_ALLHOURS_inland_event", clear


gen temp_f_drop = regexm(var,"temp_f")
drop if temp_f_drop == 1


split var, p(_)
rename var4 hour
destring hour, replace
order hour

label variable hour "Hour"

replace hour = 24 if hour == 0
sort hour 



//do it for the presentation
local color black
local style shortdash

twoway 	(scatter coef hour, color(black) lwidth(medthick) ) ///
		(line coef hour, lcolor(black) lwidth(vthin) ) ///
		(line ci_lower hour , lcolor(`color') lpattern(`style')) ///
		(line ci_upper hour , lcolor(`color') lpattern(`style')), ///
		yline(0, lcolor(black))   graphregion(color(white)) bgcolor(white) ylabel(,grid glcolor(gs15) angle(0) ) ///
		ytitle(Coefficient, height(-70)) graphregion(color(white)) bgcolor(white) ///
		xtitle(Hour, height(5)) ytitle(, placement(3) ) ylab(,angle(horizontal)) ysc(titlegap(-10) outergap(0))  ///
		legend(order(1 3) label(1 Hourly impact) label(3 95% confidence interval)) ///
		xlab(2 "2am" 5 "5am" 8 "8am" 11 "11am"  14 "2pm" 17 "6pm" 20 "8pm" 23 "11pm")  xline(14 17, lwidth(mediumthick) lcolor(gs10)) 

graph export "latex//All_hours_inland_color.pdf", replace

